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Abstract: This article presents the results of a quantitative study of the small- x 
data at HERA, using the CCFM equation. The first step consists of choosing the 
version of the CCFM equation to be used, corresponding to selecting a particular subset 
of next-to-leading-logarithmic corrections — the choice is constrained by requiring a 
phenomenologically reasonable small- x growth. For the time being, the parts of the 
splitting functions that are finite at z = have been left out. We then examine 
results for F£ , R, the transverse energy flow, the charged-particle transverse-momentum 
spectrum and the forward-jet cross section and compare to data. While some of the 
data is reproduced better than with DGLAP-based calculations, the agreement is not 
entirely satisfactory, suggesting that the approach developed here is not yet suitable 
for detailed phenomenology. We discuss why, and suggest directions for future work. 
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1. Introduction 



The HERA accelerator has opened up a radically new region of phase space for Deep 
Inelastic Scattering (DIS). This is the region where Bjorken-x is very small, while the 
photon virtuality Q 2 is still hard but moderate. In this kinematical domain it is no 
longer clear that DGLAP evolution [1], which resums logarithms of Q 2 , should lead to a 
good description of the physics. Indeed one expects that the resummation of logarithms 
of 1/x should become equally if not more important. One of the aims of the HERA 
experiments is to study observables which have different characteristics according to 
the kind of resummation performed, in order to clarify one's understanding of the 
kinematical domain and of the details of the resummation of logarithms of x. 

For total cross sections (i.e. DIS structure functions), leading logarithms of 1/x 
(al l \n n l/x) can be resummed using the BFKL equation [2]. This predicts that the 
total cross section should grow as a power of 1/x. The rapid rise of the F 2 structure 
function at small x [3,4], when first discovered, was optimistically hailed as being 
evidence for this. There are however difficulties with such a claim at a quantitative 
level: 

1. With DGLAP-based approaches the data can be quite straightforwardly fitted, 
using either an input parton distribution which rises at small- a;, or by starting 
the evolution from a low scale [5,6]. 

2. Leading logarithmic (LL) BFKL evolution leads to much too steep a rise of F 2 
to allow a fit to the data (see e.g. [7]). 

3. The recently calculated next-to-leading (NLL) corrections to the BFKL kernel 
[8,9] turn out to dominate over the leading contribution for any realistic value of 
the strong coupling (so perhaps it should come as no surprise that leading-order 
BFKL cannot fit F 2 ), and if taken literally, lead to nonsensical results [10]. 

4. BFKL evolution is not ordered in transverse momentum, as a result of which 
diffusion in transverse momentum [11-13] leads to significant dependence on the 
infra-red region, where the perturbative series is irretrievably out of control. 

Points two and three are closely connected, and in due time, perhaps following an 
approach similar to that of [14] , we will probably learn how to resum large parts of the 
NLL corrections. One can to some extent avoid the problems of the infrared region by 
including small-x resummations in the Q 2 -evolution of the F 2 structure function [15]; 
the tradeoff is that as in DGLAP, one needs an input parton distribution over the whole 
range of x, possibly leaving a little too much freedom in fitting the data. 

A way of observing BFKL phenomena while avoiding sensitivity to poorly known 
input parton distributions or to the infrared region was suggested some time ago by 
Mueller and Navelet [16] in their proposal to study more exclusive quantities, in partic- 
ular the cross section for the production, in hadron-hadron collisions, of two similarly 



hard jets separated by a large rapidity interval, which should grow as an exponential of 
that rapidity interval. The advantage of this observable (and of many other similarly 
exclusive measurements which have since then been proposed), is that by triggering on 
two large transverse momenta at the extremes of the evolution, one can suppress the 
probability of diffusion into the infra-red region, leaving a quantity which can, at least 
in principle, be calculated perturbatively (of course if the rapidity interval between the 
two transverse momenta becomes too large, then diffusion once again becomes impor- 
tant). A second advantage is that such observables can discriminate between BFKL 
and DGLAP evolution more effectively than the total cross section at small x: DGLAP 
evolution cannot possibly mimic the expected BFKL behaviour, because generally the 
two transverse scales are of the same order, and DGLAP only enhances cross sections 
for interaction between objects with very different transverse scales. 

So an efficient way to identify the most important characteristics features of lnx 
resummation is to study the final state. It turns out that the multi-parton distribution 
of an event cannot be reliably obtained from the BFKL equation. The fundamental 
additional element that must be taken into account (as in almost any calculation of 
final-state properties [17]) is that of the coherence of QCD radiation [18]. For the 
particular case of an evolution in which the transverse scale increases along the direction 
of evolution, the correct formulation was derived in [19] and is known as the CCFM 
equation — it correctly resums the leading Inl/x and InQ terms both for inclusive 
and non-inclusive quantities [20]. Angular ordering should also reduce the dependence 
on the infra-red region, since large drops in the transverse scale are disfavoured, thus 
limiting diffusion [21]. 

In contrast to the BFKL and DGLAP equations, whose parton densities depend on 
two variables (energy fraction and transverse momentum), the CCFM gluon density 
depends additionally on a limiting angle. This increases very considerably the math- 
ematical complexity of the equation, and it can only be solved numerically (and even 
that is highly non-trivial). 

In view of this difficulty one is led to ask whether it is really necessary to use the 
CCFM equation in order to get a valid description of the final state. It has been 
pointed out by Ryskin [22] that if one formulates the BFKL equation as an evolution 
in rapidity rather than in x, then one automatically has angular ordering. This for 
example is the approach that has been taken by Orr and Stirling [23] in calculating 
Mueller-Navelet-like two-jet observables at the Tevatron, and in such a context it may 
be valid. 

However the point of the CCFM equation is not simply that of angular ordering - 
rather it consists in the combination of the ordering of angles and the ordering of the 
x variables, in a situation in which transverse momenta increase along the chain, i.e. 
a DIS-type situation. If one applies the BFKL equation with only rapidity ordering 
to a DIS situation then one will obtain an answer for the cross section with spurious 
double transverse logarithms (a 2 In 1/x In 2 Q) at NLL order. This problem is closely 



connected with that of the choice of scale in the kernel [9, 14] — indeed in the language 
of [9], the BFKL equation with rapidity ordering corresponds to choosing a symmetric 
scale kik2 in the LL kernel. What this means in practice is that the cross section at 
small x and large Q 2 , instead of behaving, as it should do, as 



It turns out though, that even the BFKL equation with ordering in x (which does not 
suffer from the above problem) manages to reproduce correctly certain properties of 
the CCFM final state: in a fixed-order calculation by Forshaw and Sabio Vera [24], 
subsequently extended to all orders by Webber [25] it was demonstrated that in the 
double-logarithmic approximation all one-particle inclusive quantities have the same 
leading-logarithmic terms in both approaches. Nevertheless one should be wary of such 
a result, since there is the distinct possibility of one-particle-inclusive quantities having 
a collinear divergence at subleading order in BFKL [26], and secondly, BFKL may not 
correctly reproduce correlations between emissions [26,27]. 

So it is likely that the CCFM equation is the best approach currently available for 
studying the final state. At the same time it should be noted that the CCFM equation 
is not entirely free of problems. Its approximation of ordering simultaneously in angles 
and x is valid when the transverse scales increase along the direction of evolution. 
When the transverse scales decrease, the relevant approximation is that of ordering 
simultaneously in angles and x. One of the characteristics of small-x evolution is that 
even if, predominantly, the transverse scale increases (as in DIS), there are evolution 
steps in which it drops (diffusion) - - and the CCFM equation is subject to large 
subleading corrections in those steps (again this is related to the issue of the choice of 
scales). One might be tempted to argue that in DIS the problem arises only occasionally 
and so will not matter excessively — but several final-state "BFKL-signals" at HERA 
tend to favour evolution paths in which the evolution has no dominant direction in 
transverse momentum. 

Other problems that arise are connected with two further subleading elements of 
the CCFM equation: the treatment of factors of (1 — z) in the kinematical variables of 
the virtual corrections is not defined at LL order, and as we will see, a modification of 
their treatment can lead to very large subleading corrections. Secondly the part of the 
splitting function which is finite as z — > 0, 
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can also give important sub leading corrections. 1 There are several possible ways of 
including it in the evolution, and they may well give significantly different subleading 
corrections. Because of the enhancement at z — 1, these contributions are quite difficult 
to study numerically. It should be noted that they (in particular the 1/(1 — z) part) 
are vital for a correct description of infra-red sensitive final-state properties. 

In this uncertain theoretical context, one thing is clear however: experimentally, 
DGLAP-based approaches (whether Monte-Carlo event generators or NLO calcula- 
tions) are incapable of describing many of the small- £ final-state observables that are 
supposed to be sensitive to BFKL phenomena [28-31]. So there is an urgent need for 
some BFKL-based phenomenology, and in particular a Monte Carlo event generator 
based on a leading In x resummation. This article is in some sense intended as a study 
of the extent to which one can carry out quantitative phenomenology using the CCFM 
equation as it stands currently. 

At this point, it is worthwhile specifying some realistic aims. Quite a significant 
issue is whether one should worry about matching the NLL corrections calculated for 
the BFKL equation [8,9]. If they were small, the answer would obviously be "yes." But 
since they are so large, all that should matter in practice is the resummed corrections. 
Some idea of what a (partially) resummed result might look like for the BFKL equation 
is given in [14]. Beyond this we are limited to being guided by the experimental data. 
In particular, we aim to reproduce the total cross section as a function of x and Q 2 . 
This will force us to deal with problems associated with the treatment of the infrared 
region. 

Bearing in mind these requirements, our approach is as follows: we choose a CCFM 
evolution equation which phenomenologically seems to have a reasonable exponent 
(i.e. growth at small x). By "phenomenologically reasonable," we mean that for the 
relevant values of a s it should be roughly consistent with the growth of F 2 , and with 
expectations from a resummation of NLL terms of the BFKL kernel (as for example, at 
an approximate level, in [14]). The methods used to determine the exponent follow on 
from those of [21]. Because of the lack of a well-defined way of including the z — > finite 
part of the splitting function, and because of the considerable difficulty in implementing 
each new way that one might invent, we simply leave it aside, for possible inclusion at 
a later stage. 

We then develop a suite of programs allowing us to study structure functions and 
one-particle inclusive quantities, and constrain the treatment of the infra-red region 
by fitting to the F 2 structure function data. At this point we should be in a position 
to give relatively parameter-free predictions for other structure functions (F£, F L ) and 
various (almost) one-particle-inclusive final-state quantities: the transverse energy flow, 
charged particle spectra and the forward-jet cross section. 

Before entering into the details it should be noted that there are already in exis- 
tence a number of programs based on the CCFM equation, for studying small-x DIS 

1 Here the z — > 1 singularity is regulated by a corresponding Sudakov form factor. 



A 



scattering. 

In [32] there have been studies of the F 2 , F% and Fl structure functions using the 
form of the CCFM equation without the z — > finite part of the splitting function, and 
with all other 1 — z factors replaced by 1 (both in the real and virtual corrections) . 

The program SMALLX [33] is a direct Monte Carlo implementation of the CCFM 
equation. It includes the z — > finite part of the splitting function (unlike the results 
presented here), and additionally, the exact kinematic constraints. The main problem 
with this program is that it is a forward-evolution event generator, with weights that 
have extremely large fluctuations. This renders it somewhat unwieldy to use. Ad- 
ditionally, though it contains all the physics, there is no way of modifying the exact 
implementation of the subleading corrections — such modifications turn out to have 
large effects. The SMALLX program has recently been used and studied in [34], for 
the case of the production of charm quarks — unfortunately the precision of the data 
for charm quark production is rather poor (the tagging of open charm is not a simple 
task), and there is little final-state data available. 

Another major numerical study is based on the Linked Dipole Chain (LDC) model 
[35], which has been implemented as a Monte Carlo event generator [36]. Formally it is 
based on the CCFM equation with some simplifications and modifications. Perhaps the 
most important modification is that the evolution has been made symmetric (i.e. it is 
equally valid whether the transverse scale increases or decreases along the direction of 
evolution). A lack of symmetry is one of the main defects of the CCFM equation. On 
the other hand, the LDC model has the problem that it rearranges the real and virtual 
terms of the evolution in a manner which alters the LL small-x evolution. Whether 
this is a serious problem is to some extent a question of opinion — after all, given that 
the known NLL corrections are so large, what really matters phenomenologically is the 
resummed evolution, not its perturbative expansion. The LDC approach gives satis- 
factory results for the structure function data and some final state properties such as 
the transverse energy flow. However it has problems with the charged-particle trans- 
verse momentum spectrum and with the forward-jet cross section. As will be seen, 
our approach has rather similar problems — a discussion of this will be given in the 
conclusions. 

The structure of this paper is as follows: in section ^| we review the structure of 
the CCFM equation, discuss approximations that are commonly made when solving it, 
together with the treatment of subleading corrections in the non-Sudakov form factor. 
In section [3] we examine the asymptotic exponents of various versions of the CCFM 
equation. In section |] we consider the evolution parameters and corresponding fits to 
the F2 structure function, as well as the consequences of these different fits for F 2 C and 
R. In section [5] we then discuss the determination of one-particle-inclusive final-state 
quantities, and present results. The whole approach, and directions for future work are 
finally examined in the light of these results in section ||. 

In an appendix we give details of a novel analytic approach for implementing the 



off-shell boson-gluon fusion [37-39], which is convolved with the unintegrated gluon 
structure function to give physical cross sections. Traditional approaches generally 
require a two-dimensional numerical integration over the quark- ant iquark phase-space 
(or alternatively a treatment in terms of the Mellin transform) — the solution presented 
here carries out these integrations analytically, even in the presence of massive quarks, 
greatly simplifying the determination of physical cross sections from the unintegrated 
gluon structure function. 

2. Angular ordering 
2.1. Kinematics 

In what follows, we discuss the unintegrated gluon density of a proton (taken to be 
massless) with energy E p . This gluon density, A(x, k,p), depends on: x, the fraction 
of the proton's plus-momentum component that is carried by the gluon; k the gluon's 
transverse momentum; and p, which defines the maximum-allowed angle of all emissions 
prior to this gluon, through the relation p = 2xE p tan # max /2, where E v is the proton 
energy. The recipe for choosing p in DIS will be discussed later. 

The angular ordering condition on the last emitted gluon (with transverse momen- 
tum q t and plus-momentum fraction x(l/z — 1)) can be written as 



where for convenience we have introduced the variable q — qt/(l — z). We shall also 
use the rapidity r\ of an emitted particle: 



where, q± = go ±93 and the proton has a positive 3-component. For the purpose of the 
kinematics, emitted gluons will always be taken to be massless. 

Considering a general branching (see figure 0), in which a space-like gluon k n has 
transverse momentum k n and momentum fraction x n , and an emitted gluon q n has 
transverse momentum then the general ordering is q n+ \ > z n q n (with z n = x n /x n _i 
and q n = q t;n /(l - z n )), with k n = fc n _i - (1 - z ri )q ri . Additionally all the q t>n satisfy 
q tt n > /-£ where \i is our collinear cutoff. 

The 1 — z n factor is necessary for the exact reconstruction of the transverse momenta 
(alternatively, one can interpret it as being necessary for exact angular ordering). It 
gives a next-to-leading contribution to the structure function evolution (as does angular 
ordering itself). 
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Figure 1: Kinematics of the branching process: ki n is momentum coming from proton; k n 
is momentum entering into the photon-gluon fusion quark-box. 

2.2. The CCFM equation 

The equation satisfied by the unintegrated gluon structure function, A(x, k,p) is 
A(x,k,p) = A m {x,k,p) + 

/ — / -4 «s (ft) \k + q t \,q) A(>, k, q)Q(p - zq) , (2.2) 

Jx z J *q 

where a s = OsCa/tt and ql = (l — z)q. The non-Sudakov form factor (virtual correction) 
is 

\nA(x,k,q) = - [ - I ^a s (q t ')e(k-Tq>)e(q'-zq). (2.3) 
Jx z J iiq' 

It is usual to set T = 1, however at leading- logarithmic order it is equally valid to 
choose T — 1 — z. A motivation for choosing the latter form is that in the real emission 
term, one has a dependence onfc + (l — z)q = k + q t ; one can argue that in analogy, 
since q' is a rescaled transverse momentum, q' = q' t /(l — z), it is Q(k — q' t ) and not 
the more usual Q(k — q') which should appear in the form factor. More formally the 
choice T = 1 — z, at first sight, looks like it ought to cancel the contribution from all 
emissions with q t < k. 

The two choices will be discussed in more detail in section ^. 

2.3. The initial condition 

To ensure a sensible initial condition (one which leads to a collinear safe result after 
evolution), should correspond to the distribution for a Reggeised gluon. Relating 
to a normal, "non- Reggeised," initial condition, ^4( imt ), we have 

4 row i \ f 1 dz dA (init \z,k) . 1 , 
A^\xXp) = J z - dl J /z > A{x/zM. (2.4) 



The form that we take for ^4.( init ) (x, k) is 

(x, k) = Af*^- (1 - x) n e~ k2/k o , (2.5) 

with Af, n and k parameters to be discussed below. Note that A^\x,k,p) does not 
depend on p. The normalisation is such that 

xg(x, Q 2 ) = J d 2 kA{x,k,Q) , 

where g(x,Q 2 ) is the gluon density in the proton; Af is the result of the energy sum 
rule for the initial condition: 

Al = J d 2 kdxA^\x,k) . (2.6) 

2.4. The running coupling 

We use the one-loop running coupling, taking q t as its scale (in accord with the results 
obtained from the NLL calculation [9]). Since qt kinematically can go down to small 
values, it is necessary to introduce a prescription for the running of a s at low scales: 
we modify a s so that it reaches a plateau in the infra-red region: 

with 

llC A -2n f 
bo ~ 12^ ' 

and Ao chosen so as to ensure that a s (0) = ai°\ with oi°^ a parameter; finally, we take 
n f = 4 and A = 0.135 GeV. 



2.5. Parameters 

To summarise, the parameters that must be chosen are the following: af, the value 
of a s at which it freezes, a cutoff on emitted transverse momenta q t > (and q' t > /i 
in the virtual corrections), and the parameters which affect the properties of the initial 
condition, Af, n and k^. They will be constrained by requiring a good fit to the F2 
structure function, as discussed in section £|. One consistency check of our approach 
will be that it is possible to find a set of fit parameters that is in accord with one's 
physical expectations. For example the normalisation should be such that the energy 
sum rule should integrate to Af ~ 0.5; k is expected to be of the order of 1 GeV; ai°^ 
should be of the order of 0.5; and n ~ 4. 

Parameters that are only weakly constrained by F2 should then be varied to probe 
their effect on final state quantities. 
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Figure 2: The asymptotic exponent of the BFKL equation and various forms of the CCFM 
equation, as a function of a s . The "BFKL + NLL" curve is the resummed result (scheme 3) 
of [14]. 



3. Asymptotic exponents 

Before engaging in a fit of the F 2 structure function, it is helpful to study, as a function 
of fixed a s , the "asymptotic exponents," u, of the various forms of the CCFM equation, 
since this will indicate which one is most likely to be in accord with the phenomenology 
(though it should be borne in mind that u>(a 8 ) is not physically observable). The 
technique adopted is based on that developed in [21]: 

w = -jt— lnA(x,k,p) , (3.1) 
a In a; 

with x sufficiently small that the result is independent of x, k, and p. In the BFKL 
case, uj = a s 4 In 2. 

Figure |2] shows the asymptotic exponents of several versions of the CCFM equation, 
compared to the BFKL equation, as a function of a s . The three versions of the CCFM 
equation are as follows: the case "1 — z — > 1" corresponds to ( |2.2|) and (|2.3| ) with all 1 — z 
factors replaced by 1 (also in the relation between q and q t ) . This is an approximation 
quite commonly used (studied for example in some detail in [21]) because it greatly 
simplifies the numerical solution. However it leads to a slight violation of the angular 
ordering (or to an incorrect reconstruction of the transverse momenta, depending on 
whether one interprets q as a real, or a rescaled transverse momentum). The other two 
versions of the CCFM equation implement the correct treatment of angular-ordering (or 
transverse-momentum reconstruction), but differ in their treatments of the form factor. 



The T = 1 case corresponds to the version of the CCFM equation more commonly 
quoted. 

One sees immediately that the 1 — z — > 1 and T = 1 cases have exponents which 
are rather similar to that of the leading-logarithmic BFKL equation. In particular 
one sees that for the values of a s which are typical in HERA physics, namely 0.2-0.3, 
the exponents are far above anything observed experimentally. In a non-asymptotic 
situation (i.e. x not infinitely small) these exponents can be reduced by the introduction 
of a large collinear cutoff, but, it turns out, not sufficiently to allow a fit to the data. 

So the version of the CCFM equation that we will concentrate on in the remainder of 
this article will be that with T = 1—z. One can calculate numerically its next-to-leading 
correction to the exponent, which turns out to be extremely large: us = 4 In 2a s — (75 ± 
4)«g +0 (af). This is to be compared with uj = 41n2a s — 17.91«g in the case of the full 
NLL BFKL corrections (jif = 0) [8,9]. One is immediately induced to ask why setting 
T = 1 — z produes such a large correction. Qualitatively, it leads to large negative 
virtual corrections from q' t < k for z close to 1. However the corresponding real-emission 
contribution, q t < k for z close to 1, turns out to be dynamically suppressed because 
it corresponds to a large emission angle, which reduces the phase space available for 
subsequent emissions. Consequently the asymptotic exponent receives a large negative 
correction. 

Another question is whether it matters that the NLL correction is so much larger 
than the true one. One can argue not, since the pure NLL correction turns out to be 
visible only for extremely small values of a s , while phenomeno logically we are interested 
in the region a s ~ 0.2, where formally more subleading corrections (e.g. NNLL etc.) 
are just as important. So what matters should be the resummed asymptotic exponent, 
and from figure |2| it seems that this quantity is roughly in line with phenomenological 
expectations. In fact if one carries out a partial resummation of subleading corrections 
to the exact NLL BFKL exponent (the curve shown in figure corresponds to scheme 
3 of [14]), one obtains a result which is not so different from that with the equation 
being used here. 

Having completed the discussion of the asymptotic exponent, it is however vital to 
remember that much of the HERA data is to be found in regions of x (or energy) which 
may be quite far from asymptotic. 



4. Structure functions 

4.1. Convolution with the quark box 

To obtain a result for the F 2 structure function (or any other cross section) it is neces- 
sary to perform the convolution of the unintegrated gluon density A(x, k,p) with the 
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boson-gluon fusion matrix element (which is discussed in the appendix): 




J X 



In a similar way, one can determine Fl and _F 2 C . The maximum angle is taken, as in 
the program SMALLX [33], to be defined by the rapidity of the combined qq pair 



Modifying this limit would alter our results at next-to-leading order. The light quarks 
are taken to be massless, while for the heavy quarks we use m c = 1.5 GeV and m b = 
5.0 GeV. The convolution is performed in a frame in which the proton and photon 
are collinear and along the z-axis (this differs from the laboratory frame by a rotation 
and a boost). Correspondingly, transverse momenta for final-state quantities are also 
determined in such a frame. 

As discussed in the appendix, a is determined entirely analytically. This is in 
contrast to other approaches, which generally leave a in the form of a double integral 
which must be evaluated numerically [37-41]. The reduction in the number of numerical 
integrations simplifies considerably the implementation of the convolution. The only 
disadvantage is that the analytic integration must be performed with fixed a s — so the 
value that we use is given by the logarithmic mean of a s between scales ^Jk 2 + m 2 and 



In the limits k -C Q and k 3> Q this gives the correct overall dependence on the running 
coupling (possibly incorrect with respect to details of the treatment of the quark mass). 

4.2. Structure function results 

The fit for the structure function is performed in the kinematic range x < 1CT 2 and 
8 < Q 2 < 140 GeV 2 . The parameters to which the fit is most sensitive are /i, the 
collinear cutoff, and ao, the value at which a s freezes. We find that 0.01 < \i < 0.1 GeV, 
and «o ~ 0.5 lead to reasonable fits of the HI and ZEUS 1994 structure function 
data [3,4]; as one would expect, there is relatively little dependence on the value of n 
(the power of (1 — x) in the initial condition, see Q2.5|) ), though it does affect the fitted 
value of Af. It turns out that for the fc-dependence of the initial condition, oc e~ fc2//fc o, 
there is a correlation between the value of ko and the optimal value of fi. Good fits are 
obtained with a range of k : we will use ko = 1 GeV, and then examine the effect on 
other quantities of varying it. Modifying ko by a factor of 2 (and changing [i and ao 
accordingly) affects the normalisation by about 50%. 
The set of parameters that will be used is as follows: 
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Figure 3: Fits to the F% structure function, compared with HI [3] and ZEUS [4] 1994 data 
points. 
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A comparison with the data is shown in figure [| and is seen to be in good agreement. 
The x 2 are 64.4 for 59 HI points [3] and 84.7 for 80 ZEUS points [4]. 

Using the above fit parameters one can give expectations for the charm structure 
function and R, the ratio of longitudinal to transverse structure functions. The uncer- 
tainty on the charm structure function is gauged by varying the charm mass between 
1.3 and 1.7 GeV. The results are shown in figure |]. There is some indication that our 
curves are low compared to the data. 

The results for R, shown in figure [5], are concentrated around R ~ 0.2 fairly inde- 
pendently of x and Q 2 . This seems to be fairly typical of results from BFKL-based 
calculations (see for example [44]). 





Figure 4: Results for F 2 C compared with HI 1994 data [42], and ZEUS 1994 data [43]. 

Before moving on to a discussion of final-state quantities, it is perhaps worth com- 
menting on some rather unexpected properties of the T = 1 — z form of the CCFM 
equation: as one decreases the cutoff fi, the steepness of the small- x growth decreases. 
The reason for this unusual behaviour, is that increasing \x removes from the virtual 
corrections a region which was already dynamically suppressed in the real emissions (see 
the discussion in section ^ thus eliminating a negative contribution to the evolution, 
and increasing the small-x growth. On the other hand, as is more normal, increasing 
ai°^ increases the small- x growth (though not as substantially as in the pure BFKL 
case). As a result the best fits tend to favour a diagonal band in the cuo, hi/x space. 

This is to be contrasted with the behaviour that one sees with the BFKL equation 
or with the CCFM equation with T = 1, where the small-x growth increases as one 
decreases \x (there is none of the dynamic non-cancellation between real and virtual 
parts which is seen with T = 1 — z). We note that if one tries to fit the structure 
function using BFKL, or CCFM with T = 1 evolution, then to suppress the small-x 
growth sufficiently to fit the data, one requires an unphysically large collinear cutoff 
H ~ 3 GeV. 

5. One-particle inclusive quantities 

A number of final-state properties measured at HERA can be approximated by one- 
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Figure 5: R 



particle-inclusive quantities. The latter are relatively straightforward to calculate with 
the CCFM equation once one has in place an approach for solving for the evolution of 
the unintegrated gluon distribution. 

Suppose that one is interested in determining the number density of particles en- 
tering into a certain region of phase space with a given transverse momentum q t and 
rapidity fj. It is convenient to introduce an intermediate gluon density Bq tf j(x,k,p) 
which is obtained by considering configurations with any number of emissions, followed 
by an emission into the region of interest: 



Bqtv( x i k iP) = [ — [ ^-4«s(g<) A(x/z,\k + q t \,q) 
J x z J irq z 



x A(z, k, q) Q(p - zq) 5(q t - q t ) 5 ( 77 - In 



zq 
2xK 



(5.1) 



To obtain a full one-particle-inclusive density Q £ ^(x, one then allows any number 



1 A 



of further emissions, 



C qtT) (x, k,p) — B qtV (x, k,p) + 
f 1 dz f d 2 q _ 

/ — / — 2 a s {qt) Cg tfj (x/z,\k + q t \,q) A(z,k,q) Q{p - zq) . (5.2) 

Finally one performs a convolution with the boson-gluon fusion matrix element, as in 
( f4.1| ), to obtain the single-particle-inclusive differential cross section: 

a An°A- a- = [ — [ ^a{z,k 2 ,Q 2 ,m 2 )C m {x/z,k,p) , (5.3) 
dx dQ z dqt arj J x z J k 

where in a one sums over F 2 and 2xF%, and over quark flavours. The contribution to 
the final state from the quark box is not included. 



5.1. Transverse energy flows 

The mean transverse energy flow is given by 

dE t / da \ 1 I" ^ da (5 4) 

dr) ' \dx dQ 2 J J dx dQ 2 dq t drj 

Results are compared with HI data [28] in figure || One sees that they are uniformally 
low. There are probably two main reasons why this is so: firstly we have neglected soft 
radiation from the i-channel gluon, namely the 

dz 
1 - z 

part of the splitting function. It is responsible for the bulk of the multiplicity at 
small transverse energies. As such it is a formally subleading term. However given 
that "small transverse energies" may mean of the order 1 GeV, and that the actual 
transverse energies that one is observing are about 2 GeV, one can immediately see that 
soft radiation may contribute significantly. Another element comes from hadronisation, 
at a level analogous to a 1/Q correction in e + e~ or DIS event shapes, from which one 
might expect a contribution of the order of 0.5 GeV per unit rapidity (though this 
amount is tightly correlated with what one takes as the perturbative contribution [45]). 

In association with these ideas, it is interesting to note that if one simply adds 
1 GeV of transverse energy per unit rapidity to the curves obtained with ( |5.4j ) then one 
finds a somewhat better (but not perfect) agreement with the data. The agreement 
breaks down in the highest rapidity bins at higher Q 2 values — this is to be expected 
given that we are not including the transverse energy that comes from the quark box. 

Using other parameter sets that are consistent with the structure function data has 
little effect on the E t flow. 
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Figure 6: E% in various bins of x and Q 2 ; A comparison is made with HI 1994 preliminary 
data [28]; the pseudorapidity is calculated in the photon-proton centre-of-mass system, with 
the proton direction being at negative rapidities. 



5.2. Particle spectra 

To avoid the problems associated with the E t flow, one should consider a measurement 
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Figure 7: The charged hadron fcf = q^ h spectrum in various bins of x and Q 2 for 0.5 < 
Vcms < 1-5. The data are from the HI collaboration [29]. 



which concentrates on high-transverse momenta, and is consequently less sensitive to 
hadronisation and to the mistreatment of the emission of soft particles. Such an observ- 
able is the charged hadron transverse momentum (qt,h) spectrum for a given rapidity. 



Since the one-particle inclusive cross section calculated in (|5.3|) is for emitted gluons, 
to allow a comparison with the data it is necessary to perform the convolution with the 
appropriate fragmentation function to obtain the spectrum for hadrons. The spectrum 
for a hadron of type h is 

dn u q2\ = ( d(J V 1 f d ^ D h (g tn /g t ,f) da (5.5) 
dq t ,hdr] ' \dxdQ 2 ) J qt h q t 9 *' * dxdQ 2 dq t drj 

We make the approximation that the direction of the resulting hadron is the same as 
that of the parent gluon and use the NLO fragmentation functions of Binnewies et al, 
as given in [46]. At our level of accuracy, it would have been equally valid to use the 
LO fragmentation functions. 

We compare with the preliminary results from the HI 1994 data [29] in figure [?J 
At low transverse momenta the spectra that we obtain are much lower than the data 
- this is as expected given the lack of soft-radiation and hadronisation corrections 
discussed in the previous section. 

For the large-g^ part of the spectrum there is moderate agreement for a range of x 
and Q 2 bins, though the CCFM results seem systematically low. At larger x this can be 
explained by the fact that we are not counting the charged particles that come from the 
fragmentation of the quark box. At lower x one can however say that the CCFM results 
are favoured compared to those from a DGLAP-based Monte Carlo event generator, 
such as LEPTO [47], though they are still not perfect. 

The sensitivity of these results to the evolution parameters is moderate: altering ko 
by a factor of two (and the other fit parameters appropriately) affects the normalisation 
by about 25%. 

5.3. Forward jets 

Another quantity that is supposed to be particularly sensitive to small-x dynamics is the 
forward-jet cross section, originally proposed in a slightly different form by Mueller and 
Navelet for proton-proton collisions [16] and then adapted to DIS [48,49]. One of the 
main difficulties in obtaining a reliable prediction of this quantity is that experimentally 
one measures a jet, using a jet algorithm. The best that one can do with a one-particle- 
inclusive cross section is to associate with an emitted gluon a jet with the same energy 
and direction. This neglects the effect of multiple emission from the t-channel gluon, 
partonic showering and hadronisation, all of which play a role which varies according to 
the jet algorithm: as shown by the ZEUS collaboration [30], the choice of jet-algorithm 
can quite easily affect the results by about 15%. So in some sense this is the highest 
level of precision with which one can compare to a single-particle-inclusive approach 
with the data. Another element is that events in which 2 jets satisfy the selection cuts 
are treated differently experimentally (one counts only one of them), and one-particle- 
inclusively (one counts both of them) - however such events occur rarely (2% of the 
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time in the ZEUS sample), and so this difference in treatment is of little practical 
importance. 

In most forward-jet rate calculations, rather than using CCFM (or BFKL) evolution 
before the emission of the forward jet, as in (|5.1[) , one uses DGLAP evolution (i.e. rather 
than B depending on A, it depends on the normal DGLAP structure functions, xg(x) 
and xq(x)). The motivations for using DGLAP evolution are that at relevant values of 
Xj et (which can be of the order of 0.1) one has a good knowledge of the DGLAP parton 
densities; our CCFM-based densities in that region may not be sufficiently constrained 
by the small-x F 2 data, and in any case do not include quarks. On the other hand, 
with CCFM evolution one has the full (k,p) dependence of the gluon density, whereas 
with DGLAP evolution this information must be added in. 

To use the normal DGLAP structure functions before the emission, one follows 
exactly the same procedure as before, except that one replaces ( |5.1|) with 

' x dz f d 2 q 



J x z J z \ z / 



x A(z, k, q) 0(p - zq) 8(q t -q t ) 5 I 77 - In — — ) , (5.6) 



2xE, 



with q t = (1 — z)q and 

F(x,q 2 ) 



9 (x, qt) + 7^ (q(x, q 2 t ) + q(x, q 2 )) 



(5.7) 



where g and q are usual DGLAP parton densities in the proton (and the sum over 
flavours is implicit for the quark and antiquark densities). Fixing k through 8 2 {k + 
q t ) corresponds to the approximation of strong transverse- momentum ordering (which 
should be valid in the region of Xj et of interest). 

The cuts used for the forward-jet measurement tend to be considerably more com- 
plicated than for the energy-flow or particle spectra. As a result the simplest way to 
obtain a prediction is to use a Monte Carlo integration method which generates x, Q 2 
values, a four vector for the forward gluon (or quark), and a corresponding weight. To 
these one applies the experimental cuts as if the forward parton were a jet, and finally 
adds the weights of the configurations that survive. 

We will compare our calculations with the ZEUS measurement [30]. The cuts are: 
0.00045 < x < 0.045, E tJet > 5 GeV, 0.5 < E^JQ 2 < 2, the energy of the scattered 
electron E e i > 10 GeV, ysj > 0.1, 9j et < 8.5 degrees, xj et > 0.036, and the jet must be 
in the target hemisphere of the Breit frame. 

Figure |8| shows the ZEUS results and a calculation which includes CCFM evolution 
both before and after the forward-jet emission. The CCFM cross section is consistently 
too low everywhere. In the higher-x bins this is not unexpected since there are other 
mechanisms of importance there which are not being considered (e.g. the forward jet 
can come from one of the legs of the quark box). But at low x it is more difficult to 
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Figure 8: The forward-jet cross section, and the results of the ZEUS measurement [30]. 
The upper figure shows three different approaches to calculating the results using the CCFM 
equation, while the lower figure compares the CCFM results with other theoretical predictions 
[40,50]. 



justify. To examine the degree of uncertainty on our result, we have tried to simulate 
the effect of hadronisation and "mis-reconstruction" from a jet algorithm by adding 
1 GeV of transverse energy to all jets, while maintaining their rapidity. This has a 
greater effect at small x than at larger x, due to the fact that the cut on E 2 - et /Q 2 , 
in conjunction with the correlation between x and Q 2 (the cut on i/Bj) causes jets in 
events with larger x to have a larger transverse momentum and hence be less sensitive 
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to hadronisation. The effect of this is not inconsistent with the O (15%) jet-algorithm 
dependence mentioned above. In any case it is not sufficient to bring the CCFM curves 
into agreement with the data. The use of the DGLAP distribution for the incoming 
partons, rather than the CCFM distribution, tends to lower the results — so that 
does not help either. Changing the parameter sets for the evolution also had little 
effect, affecting only moderately (~ 20%) the normalisation. We note that a similar 
calculation carried out with the cuts used by HI also falls below their recent data [31]. 

For reference we include also a comparison of our results with those of [40], where 
a normal BFKL-based calculation of the forward-jet cross section was carried out: the 
data are considerably closer to the results of [40] than to ours. We have also included 
the results of a NLO calculation using MEPJET [50], as presented in the ZEUS forward- 
jet paper [30]. For this kind of observable NLO is the first order which contributes, 
so the fact that the CCFM-based results are above the MEPJET prediction is to be 
expected. 

Why does the version of the CCFM equation that is being used here lead to forward- 
jet results which are not satisfactory? The forward-jet measurement that is used ex- 
perimentally does not cast much light on this, in part because of the numerous cuts, 
which lead to correlations between x and E t j e t and restrict the available x range. So 
it is instructive to examine a forward-jet-like (where the forward-jet cannot come from 
the quark box) measurement which is free of these problems, such as 

dF 2 (x,Q 2 ) 

as a function of x, for fixed E t j e t, Q 2 and i]. At very small x we know that in the absence 
of small- x evolution (the Born approximation) it should go to a constant, independent 
of x. So what we can do is compare the CCFM result for this quantity with a similar 
quantity worked out in the Born approximation, namely using 

1 dz f d 2 q _ . x (X 



C^ RN \x,k,p)= [ *L I ^a s (q t ) -F(*,A8*(k + q t ) 



X e(p - zq) 8{q t - q t ) 5 (r} - In T^Tj ■ ( 5 - 9 ) 

The result is shown in figure 

Concentrating first on the Born result, we see that it takes a considerable range in 
x to reach its constant value (bearing in mind that Xj et — 0.06): assuming that the 
incoming gluon has transverse momentum Q, the kinematics of the quark box require 
the momentum fraction of the incoming gluon, x g , to satisfy x g > 2xsj — to obtain 
the full asymptotic cross section x g should be larger than that limit by an order of 
magnitude (the cross section goes as dx g /x 2 g ). But x g itself must satisfy x g <C Xj et 
otherwise the cross section is sensitive to the DGLAP parton distributions at x values 
larger than Xj et , and the gluon distribution drops quite fast with increasing x. 




Figure 9: A "toy" forward jet measurement, comparing the result with evolution to that 
without evolution (Born). In both cases the input parton distribution is the DGLAP one. 

What is surprising is that the full CCFM result should drop below the Born result 
for such a considerable range in x, and only for x < 10 -5 does it start to overtake 
the Born result: the CCFM evolution takes a long time to get started. This again is 
probably a consequence of the large virtual corrections, which can be compensated only 
when there is a sufficiently asymptotic gluon transverse-momentum distribution — and 
that which arises immediately after the forward-jet emission is far from asymptotic. 

6. Conclusions 

The main aim of this paper has been to apply the angular ordered small-x formalism, 
in the form of the CCFM equation, to the description of HERA data, with the purpose 
of determining whether, as it stands currently, it is suitable as a basis for quantitative 
phenomenology. 

As a first step it was necessary to decide on the treatment of 1 — z kinematical factors 
in the virtual corrections (T = 1 or T = 1 — z, in ( |2.3|) ). The final choice was made 
on the basis that T = 1 — z gives an asymptotic exponent which is considerably lower 
than T = 1 and LL BFKL, and roughly in accord with phenomenological expectations 
and with first estimates from a partial resummation of NLL corrections to the BFKL 
equation. This choice does however have the property that its NLL corrections are four 
times larger than the true NLL corrections to the BFKL equation. The part of the 
splitting function which is finite as z — > was left out, partly because it is not clear 
how it should be included, partly because whichever method of inclusion one chooses, 
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there are considerable technical difficulties in its implementation. We will return to 
the discussion of these points shortly. To constrain the infra-red properties of the 
evolution, we used the F2 structure function; we then examined some of the classic 
"BFKL" structure-function and final-state signals: F 2 C , R, the transverse energy flow, 
the charged hadron transverse momentum spectrum and the forward-jet cross section. 

The transverse energy flow is much lower than the data, but this was to be expected: 
both the z — > finite part of the splitting function and non-perturbative contributions 
are expected to give large (though formally subleading) contributions. 

The failure to reproduce the forward-jet cross section data is less acceptable. An 
analysis of the problem indicates that after the emission of the forward jet, the small- a; 
evolution has a long period in which there is no growth of the cross section — non- 
asymptotic effects are important. The exact reason for this may well be connected with 
the non-trivial interplay between real and virtual contributions, which leads also to the 
large NLL corrections to the evolution. The importance of non-asymptotic effects may 
also be responsible for the lowness of the predictions for F% and the charged-particle 
transverse momentum spectra (though the results for the latter are better than those 
from DGLAP-based approaches, at least for some values of x and Q 2 ). 

It is interesting to note that the LDC model [35, 36] displays problems similar to 
those found here. The physics contained is somewhat different (e.g. the evolution is 
symmetric) but the general approach is somewhat similar — they too use a fit to the F 2 
structure function as a constraint on the evolution parameters. The fit to the structure 
function data works well, but again a number of final state quantities are too low. It 
is not clear whether the problems in the LDC model are of the same nature as those 
in our CCFM implementation (i.e. due to non- asymptotic effects), but it might be a 
matter worthy of further investigation. 

As a consequence of the problems described here, one is essentially led to rule out 
the CCFM equation with T = 1 — z and without the z — > finite part of the splitting 
function. It is tempting to suggest that we should have used T = 1 and then included 
the z — > finite part of the splitting function (in some arbitrary way). This would 
be the physics of SMALLX program — which so far has only been used for the study 
of final states in which the photon-gluon fusion produces charm quarks [34]. Since 
most final state studies at HERA do not select for the presence of open charm, there 
is relatively little data for comparison, so it would be worthwhile extending that study 
to the case with the full mix of quarks in the photon-gluon fusion. 

In the longer term we need to work towards a theoretical approach free of the 
current ambiguities. Much of the basic information is probably already at our disposal, 
in the form of the NLL corrections to the BFKL kernel [8,9]. We should be looking 
for an evolution equation which rather just than giving a "reasonable exponent" (a 
condition satisfied both by the LDC and by the form of the CCFM equation described 
in this paper), explicitly reproduces the structure of the NLL corrections (whether 
the part associated with symmetrisation, of with the z — > finite part of the splitting 
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function), and which resums these corrections in a sensible way. Of course it should also 
satisfy basic properties associated with coherence in the final state. Another important 
element lies in the calculation of the boson-gluon fusion hard matrix element to next- 
to-leading order. After the completion of such a programme we should be in a relatively 
strong position to carry out quantitatively meaningful BFKL phenomenology. It may 
take a while for such a stage to be reached — but in the mean time it is of paramount 
importance that experimental measurements of BFKL signals continue. 
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A. Photon- gluon fusion 

In this appendix we outline the calculation of the one-loop photon-gluon fusion hard 
matrix element. The results presented have been derived with the aid of the HIP 
package for Mathematica [51]. 

The starting point is the tree-level amplitude, coming from the diagrams illustrated 
in figure |lO|a: 
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t — m 

l a [(Pi ~ k) ■ 7 + m] 7 M 



7 M [(Pi ~ g) • 7 + m] 7 a | ( ^ :| , 



u — m 2 



V(P2) ■ 



Here cq is the fractional quark charge, g s is the strong coupling constant, t a are the 
SU(3) colour generators, m is the quark mass, p\ and p 2 are the outgoing quark and 
anti-quark momenta and q and k denote respectively the photon and gluon momenta 2 . 
The Mandelstam variables are defined as s = (k + q) 2 , t = (q — px) and u = (q — P2) 2 ■ 
The hard matrix tensor H^* ,aa * is obtained by squaring the amplitude, averaging 
over the colour of the incoming gluon, summing over the final state (i.e. quark and anti- 
quark) degrees of freedom and integrating over the two-particle phase space d<p {pi,P2)'- 

H^*> aa \q,k) = j d(f)(p 1 ,p 2 ) J2 M ^ M ^ a * ■ ( A ' 2 ) 
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It depends only on the photon and gluon Lorentz indices and momenta (figure [10]). 



2 In this appendix k is the gluon four-momentum, and kt the transverse momentum, whereas in the 
rest of the paper we denote by k the latter quantity. 
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a) b) 

Figure 10: a) The one loop diagrams contributing to photon gluon fusion vertex, b) The 
squared amplitude Lorentz indices. 



The "reduced cross section" o> 2 appearing in (|4.1|) is then defined through 



2x f 



BFKL,aa* i 



(A.3) 



where the photon Lorentz indices are contracted over the usual projector used in DIS 
to extract the structure function F 2 and the longitudinal structure function F L : 
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and the gluon Lorentz indices are contracted over the projector 
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coming from the high energy (or k t ) factorisation prescription [37-39]. Additionally 
the prescription specifies that in ( |A.3[) , the incoming gluon momentum should be ap- 
proximated by k = Xfp + kj_ with p the hadron momentum, so that k ■ p = and 
h 2 — h 2 — -h 2 

After the contractions with the two projectors, the reduced cross section acquires 
a dependence on the proton momentum (i.e. it is not just a function of the photon 
and gluon momenta), which complicates the phase-space integration for the outgoing 
quarks. As a result it seems very difficult to carry out this integration analytically, 
and the answer is usually left in the form of a Mellin transform, or of a double integral 
which must be evaluated numerically [37-41]. 

In this appendix we show that it is possible to carry out the phase space integration 
analytically, but this must be done before the contraction with the projectors. Accord- 
ingly, we expand H^^*' aa * (q, k) over a tensor basis {T/^ ' aa (q, k)}, taking into account 
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the symmetry of H^*' 00 * with respect to both the /x, /i* and a, a* pair of indices, and 
its transversality with respect to q and k (i.e. q^H^*^* = 0, k a H^*' aa * = 0). 

First of all, let us introduce a four-dimensional basis hf u (a, b) for a two- index sym- 
metric tensor which depends on two four-momenta a and b: 
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(A.7) 



We then construct a basis ' aa (q, k) for a four-index tensor which depends on the 
two momenta q and k and which is symmetric with respect to /i, fi* and to a, a*. There 
are 21 such independent structures. Sixteen can be chosen as the products 



B^^j («» k ) = h f (?> k ) h T* (*. , <, J = 1, • • • , 4 , (A.8) 



while for the others we take 



n 17 




(?,*) 


-°18 


,aa* 


(?,*) 


n 19 




(9,*) 


RW* 

^20 


,aa* 


(?,*) 


-°21 


,aa* 


(?,*) 



5 



^,2 1 y 1 y 1 y / > 

_L (g^ q ^k a * + g^q^k a + g»* a q»k a * + g»* a * q»k a ) , 

(g^k^q*** + g^k»*q a + g^Wq 01 ' + 0"* a Vg Q ) , 
q ■ k ' ' 



Since H^* ,aa * (q,k) is explicitly transverse with respect to both g and fc, let us now 



concentrate on the set of symmetric tensors ,aa = Y^=i CijB^ ' aa that satisfy 



q ^ T w*, aa * = and kaT w*,aa* = g There are g guch independent tensors, specified by 
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With P=W$ 



This choice has the advantage that the corresponding "projectors" {T i:lipj * :Cta * (q, k)} : 



which satisfy the relations >aa * T w ,aa = 5^, have straightforward forms: 



if w = 1 + g2w * ) (<r + f, ' 2<r<r 



(k ■ q) 2 — k 2 q 2 J \ (k ■ q) 2 — k 2 q 



(k ■ q) 2 — k 2 q 2 J \ (k ■ q) 2 — k 2 q 2 J ' 

7>" - 1 /V + \ (V + :U ''' 2>l " >r ) (All) 

" 4 V + (k-q) 2 -k 2 q 2 ) V (A;-g) 2 -fcV j ' 



k 2 n 2 

h ~ 2{k . q)2 [9 9 ±9 9 ) ■ 



We then write the decomposition 



H w;aa* = H{T^' aa * , (A. 12) 



1=1 



where the coefficients if, are 



H l = T h ^ yaa , ' aa . (A. 13) 



Since the ifj depend only on k and g, there exists a frame, namely the photon-gluon 
centre of mass frame, in which the quark- ant iquark phase-space integration can be 
performed analytically without difficulty. This would not have been the case if we 
had already contracted with the photon and gluon projectors, as these would have 
introduced a non-trivial dependence on an additional external momentum (that of the 
proton). Our results for the coefficients Hi are: 
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a s e 2 Q 



Am 2 . J 



2 t2 



+ H„ 
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(A. 14) 



no 



where J = a/ (k ■ q) 2 — k 2 q 2 and 
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(A.15) 



(b) 



H {b) - 

n 2 — 



(b) 



H. 



(b) 



H. 



n A — 







Am 2 
s 


J 2 






Am 2 
s 


J 2 






Am 2 
s 


J 2 






Am 2 
s 



J 2 



k 2 q 2 s 2 
AJ 2 



t2 2 

J — m s — 



k 2 q 2 s 2 
AJ 2 



t2 2,2 

J — m s + q s 



J 2 + k 2 s — m z s — 



2 k 2 q 2 s 2 



AJ 2 



3k 2 q 2 s 2 
AJ 2 



J 2 + k 2 s - m 2 s + q 2 s + 



*rn? k 2 q 2 s 2 



J 4 



_ *£k 2 q 2 (2 J 2 - Am 2 s + s 2 ) 



J 2 (k- q y 



(A.16) 



on 



and 



2J 6 



8J 4 (k • q) 2 - 8J 4 (k ■ q)s - 8J 2 (k ■ q)k 2 q 2 s + 4JV + 



+ QJ 2 k 2 q 2 s 2 + 3k 4 q 4 s 2 + 8J 2 m 2 (-2J 2 m 2 + (k ■ q) 2 s) 
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2J 6 
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(A.17) 



+ UJ 2 k 2 q 2 s 2 + 9k 4 q 4 s 2 + 8J 2 m 2 (-2J 2 m 2 + 2J 2 q 2 - J 2 s + 2(k ■ q) 2 s) 
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2J 6 
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((k-q) -m 2 ), 
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Though individually, each of the H suffer from a singularity at J = 0, the results for 
the reduced cross sections are finite at J = 0, and are given by: 



a Fi (z,klQ 2 ,m 2 ) = ^- z k 2 a t [q 2 (\ - 1) - k 2 , -k 2 , -Q 2 ,m 2 ^ , 



(A.18) 
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where 



and 



-y 4m 2 



- 32J 6 - 72J 4 k 2 q 2 + 96J 4 (k ■ q)s + 



— ^ - O^U £iO nj y T ^ yj' 5 -T 

+360J 2 (£; • g)£;y s - 48J 4 s 2 - 330 J 2 k 2 q 2 s 2 - 315A; 4 g 4 s 2 



1 - 



4m 2 



2J 4 



16 J 6 + 



24J 4 k 2 q 2 - 32J 4 (k ■ q)s - lQJ 4 m 2 s - 48J 2 (k ■ q)k 2 q 2 s + 
6 JV + 3QJ 2 k 2 q 2 s 2 + 15k 4 q 4 s 2 ] , (A.20) 



-36J 2 £; 2 my s + 16 J 4 s 2 + 36J 2 k 2 q 2 s 2 + 15k 4 q 4 s 2 
= — 32J 8 + 80J 6 k 2 q 2 + 72J 4 k 4 q 4 - 32J 6 (k ■ q)s - 336J 4 (k ■ q)k 2 q 2 s + 
-360J 2 (A; • q)k 4 q 4 s + 16JV + 240J 4 A; 2 g 2 s 2 + 540J 2 k 4 q 4 s 2 + 315JfeVs 



IF 



r 2 {k ■ q)k 4 q 4 s + 16JV + 240JW s 2 + 540J 2 fc 4 gV + 315A; 6 gV 

- 16J 4 (A; • q) - 8J 4 m 2 - 36J 2 (A; • q)k 2 q 2 - 18J 2 A; 2 m 2 g 2 + 12J 4 s + 
+60 J 2 k 2 q 2 s + Ahk 4 q 4 s] , 



hf = 



1 - 



4m 2 



- 2U 4 k 2 q 2 + 16J 4 (fc • q)s + 



J 4 

120J 2 (A; • q)k 2 q 2 s - 8J 4 s 2 - 110J 2 k 2 q 2 s 2 - W5k 4 q 4 s 

I _ 4m 2 y2 

[8 J 4 q 2 + 8 J 4 s + 8J 2 k 2 q 2 s - l2J 2 m 2 q 2 s + 8 J 2 q 4 s + 



+8 J 2 q 2 s 2 + 5k 2 q 4 s 2 



2k 2 q 2 



h L = — jr 



16 J 6 + 



2AJ 4 k 2 q 2 - 9QJ 4 (k ■ q)s - 120J 2 (k ■ q)k 2 q 2 s + 

_ _o - non a An 



(A.21) 



+72JV + 180J 2 A; 2 g 2 s 2 + 105fcVs 2 



{d) 16m 2 g 2 



4J 4 - 12J 2 (A; • q)k 2 - QJ 2 k 2 m 2 + 



+18J 2 k 2 s + 15k 4 q 2 s . 
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